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TIME  DEVELOPMENT  OF  CERENKOV  RADIATION 

Fred  R.  Buskirk 
John  R.  Neiahbours 


Physics  Department 
Naval  Postgraduate  School 
Monterey,  California  93943 


ABSTRACT 

Most  developments  of  Cerenkov  Radiation  are  in  terms  of  the 
Fourier  components  of  the  fields  and  power  emitted  by  a  single 
electron.   When  many  electrons  in  a  compact  bunch  are  emitted  from 
an  accelerator,  the  bunch  radiates  coherently  and  at  a  lower 
frequency  than  for  a  single  electron.   The  theory  for  the  time 
structure  of  the  fields  arising  from  a  charge  bunch  is  developed, 
and  it  is  shown  that  the  source  of  the  radiation  is  di/dt. 
Present  detector  technology  should  be  able  to  resolve  these 
fields . 


TIME  DEVELOPMENT  OF  CERENKQV  I  \DIATION 

INTRODUCTION 

Cerenkov  radiation,  produced  by  a  charqe  or  qroup  of  charqes, 
novinq  faster  than  the  speed  of  electromaqnet ic  radiation  in  a 

nedium,  has  been  invest iqated ,  startinq  with  the  exDeriments  of 

1  .  2  . 

Zerenkov   in  1934  and  the  exDlanation  by  Frank  and  Tamm   in  1937. 

Since  power  radiated  by  a  sinqle  charaed  particle  is  orooortional 

to  the  frequency,  most  of  the  research  effort  has  been  devoted  to 

the  relatively  intense  optical  radiation  which  is  favored  over  the 

34 
nicrowave  reqion  by  a  factor  of  about  10  .   The  ootical  results  ' 

are  qiven  in  terms  of  the  Fourier  components  of  the  fields  and  the 

radiated  power. 

In  our  previous  work  '   it  was  noted  that  microwave  radiation 

:an  be  siqnificant  because  all  the  electrons  in  an  accelerator 

Dunch  (about  10 9)  radiate  coherently;   an  effect  which  more  than 

offsets  the  sinqle  particle  increase  in  radiated  cower  with 

frequency.   For  an  electron  beam  generated  by  a  travelinq  wave 

Ijinac  and  passing  throuqh  air,  it  was  shown  that  the  various 

larmonics  of  the  basic  frequency  up  to  about  the  tenth  are  emitted 

[In  the  case  of  an  L  or  S  band  Linac,  these  corresDond  to  10  GHz 

4 

ind  30  GHz  respectively.) 

The  time  structure  of  Cerenkov  radiation  fields  in  the  ootical 
ind  even  in  the  hiqher  frequency  microwave  reqions  is  difficult  to 
observe  because  the  detectors  reaister  power.   One  of  the  few 


7 

treatments  of  the  time  dependence,  by  Tamm   in  1939,  showed  that 

the  optical  radiation  by  an  electron  is  singular  on  the  Cerenkov 
front.   Here  we  consider  the  time  structure  of  fields  generated 
when  electron  bunches  radiate  coherently;  in  a  development  which 
complements  the  frequency  domain  analysis  of  our  earlier  work   ' 

The  fields  should  be  observable  for  beams  from  induction 
accelerators  which  produce  bunches  much  longer  than  those  produced 
by  S  or  L  band  Linacs. 


MAGNETIC  RADIATION  FIELD 

The  purpose  of  this  paoer  is  to  oresent  a  development  of  the 
time  dependence  of  the  electric  field  generated  by  the  Cerenkov- 
mechanism.   The  method  is  first  to  determine  the  potentials  from 
the  moving  charge  distribution,  and  subsequently  to  obtain  the 
fields  (in  cgs  units)  from  the  potentials  by 

B  =  V  x  A  (1) 


-*• 


E  "   v*  "  c~  TE 

0 


We  assume  a  charge  density  function  p   and  a  current  density 
j   =   Prv/c   with  the  velocity  v  in  the  plus  z  direction.   The 
charge  and  current  are  assumed  to  be  concentrated  along  the  z  axis 
such  that 

Pv(r,  t)  =  p(z,t)  6(x)  5(y)  (3) 


and  the  charge  is  assumed  to  move  with  no  change  in  shape  so  that 
the  z  and  t  dependence  of  the  charge  is 


p(z,t)  =  p  (z-vt)  (4) 

o 


Note  that  p   and  i   represent  the  usual  charge  and  current 
v       v 

densities,  while  p  and  p   throughout  this  paper  are  charge  per 
unit  length.   The  velocity  of  light  is  c  and  c   in  the  medium  and 


free  space,  respectively. 


The  potentials  are  found  by  taking  the  usual  retarded 
solutions  to  the  wave  equations;  which  become  under  the 
assumption  of  a  line  distribution  of  charge  (3), 

*(r,  t)  =  c~lj      R-Mr'  ,  t'  )  dz' 


(5) 


A(r,t)  =  |-  /  R-lp(r',f)  dz'  (6) 

o 


+    -»■-»• 

where  R  =  r  -  r   and  t   is  the  retarded  time 


t'  =  t  -   |r  -  r' l/c 

Now  (4),  the  assumption  of  rigid  motion  of  the  charge 
distribution,  can  be  incorporated  into  the  potentials,  and  a  new 
variable  u(z')  =  z'  -  vt '  can  be  introduced  so  that  the  potential 
(5)  and  (6)  become 


Mr,  t)  =  e~l    /  R_1p  (u)  dz'  (8) 


A  (r,t)  =  ^   /  R-!p  (u)  dzl  O) 

c  '  o 

0 


Also,  since  the  charge  is  confined  to  the  z'  axis,  the  new 
variable  u(z')  can  be  written  more  explicitly. 


u(z')  =  z'  -  vt  +  -  [x2  +  y2  +  (z  -  z')2] 

c 


1/2       (10) 


The  magnetic  field  B  may  be  calculated  from  (1)  and  since  A 

has  only  a  z-component,  B  has  only  the  x  and  y  components, 

3  3 

B   =  -r—  A   and  B   =  -  — -  A  .   Carrying  out  the  differentiation 
x    3y   z       y      3x   z        2       ^ 

for  the  x  component  gives 


v       3    -1 


B   =  —   /   ■£-  R    p  (u)dz' 
x    c    J       3y       o 
o 


(id 


+  V~    I  R"'  17  po(u)  dz' 

0 

For  radiation,  the  first  integral,  falling  off  as  R~ 2  at  large 
distances,  will  be  neglected  and  only  the  second  term  will  be 
considered  further.   From  (10),  it  is  seen  that  u  is  a  function 
of  x  and  y  so  that  the  second  integral  can  be  written 

v  2       y 

B   =  /   —   p  '  (u)dz'  (12  ) 

cc       Rz 
o 

where  p  '(u)  is  the  derivative  of  p   with  respect  to  its 
o  o 

argument  u.   The  corresponding  expression  for  B   has  y  replaced 
by  (-x).   These  two  components  can  be  combined  to  give  the 
total  magnetic  radiation  field  B.   In  the  cylindrical 

coordinates,  (s,  9,  z)  where  s  is  the  radius  vector 

1  /2   + 
s  =  (x2  +  y2)    ,  B  is  tangential  (i.e.  in  the  8  direction) 

with  a  magnitude  given  by 

v 2     s 


B  =  J  —   p  '  (u)dz'  (13) 

o 


cc     R2    ° 


TI^E  DEVELOPMENT 

In  order  to  evaluate  (13)  for  B,  it  is  necessary  to  consider 

the  dependence  of  u  on  z"  as  qiven  in  (10).   in  the  u-z'  plane, 

the  first  two  terms  are  a  straiqht  line  with  unit  slope  and  an 

intercept  which  chanaes  with  time,  while  the  third  term  is  a 

hyperoola  openina  in  the  +u  direction  with  asvmPtotic  slopes  of 

±  — .   The  sum  of  these  two  curves  is  u(z').   In  the  Cerenkov 
c 

case  with  v  >  c,  the  result  is  a  curve  whose  ends  both  point 
upward  as  shown  in  Fiq.  1.   As  time  increases,  the  entire  curve 
will  translate  downward  to  smaller  u  values  as  a  result  of  the 
neqative  second  term  in  (10). 

Only  chanqinq  currents  (those  with  a  non  zero  p  '  will 

o 

contribute  to  the  maqnetic  radiation  field  (13).   To  proceed 

and  demonstrate  the  method,  a  ramo-front  current  pulse  is 

chosen  as  a  simple  example.   Assuminq  that  the  front  end  of  a 

current  pulse  increases  linearly  up  to  a  constant  value,  the 

derivative  p  '(u)  will  be  a  constant  valued  square  pulse  of 
o 

maqnitude  p  '  as  is  also  shown  in  Fiq.  1.   The  correspondinq 
m 

neqative  p  '(u)  pulse  occurrinq  at  the  tail  of  a  current  pulse 

is  not  shown  and  its  effect  is  considered  separately. 

For  large  neqative  times,  the  u(z')  curve  (a)  is  completely 

above  the  oulse-like  non  zero  portion  of  the  p  '(u)  curve  so 

o 

that  the  contribution  in  (13)  to  B  from  p  '(u)  is  zero  and 

o 

therefore,  B  is  zero.   As  time  increases,  the  u(z')  curve 
moves  downward  until  the  B  oulse  begins  when  u(z*)  is  tanqent 
(curve  b)  to  the  upper  portion  of  the  p  '(u)  pulse.   The  value 


of  the  integral  in  (13)  increases  as  u(z)  continues  its 
constant  downward  motion  with  increasing  time  until  u(z'  ) 
becomes  tangent  (curve  3)  with  the  lower  part  of  the  p  ' (u) 

0 

pulse.   At  this  time  the  non  zero  part  of  the  integral  has  the 

largest  extent  -from  Zj  to  z2«   At  later  times,  the  integral 

breaks  into  two  regions  of  the  z'  axis  and  if  p  ' (u)  is 

constant,  the  value  of  the  integral  decreases  with  increasing 

time  because  the  extent  of  the  integral  in  the  two  regions 

continues  to  decrease  as  a  result  of  the  upward  turn  of  u(z'). 

This  calculation  may  be  carried  further  to  determine  the 

time  structure  (shape)  of  the  resulting  B  pulse.   Although  the 

expression  (13)  for  B  can  be  integrated  directly  in  the  case 

where  the  slope  p  ' (u)  is  constant,  it  is  instructive  to  carry 

o 

out  the  calculation  by  developing  u  in  a  power  series. 

Denoting  z  '  as  the  value  of  z'  at  which  u(z')  has  zero  slope, 
o 

the  values  of    z  ',  u(z  ')  and  the  second  derivative  are 

o       o 


v2         »   -1/2 


z  J  =  -  s   [  —    -  1)  (14) 


c 


V2      \ -1/2 


u(z  _')  =  s   —  -  M      -vt  (15) 


92U     I       1  C2  /V2      v   3/2 
=-  — 1  =2A        (16) 

3z*2   i  z     s  v2  V  c2     / 
o 


so  that  u  can  be  expressed  as  a  power  series  about  the  minimum 


u  =  u(z  ' )  +  A(z'  -  z  ' ) 2  (17) 

o  o 


The  limits  z ,'  and  z 2'  can  be  written  in  terms  of  the 

minimum  value  as  z ' ?  =  z  '  +  Az  and  zi'  =  z  '  -  Az  '  where  Az  ' 

*     o  -1      o 

is  the  value  of  z ' -z  '  such  that  the  difference 

o 

u(z')  -  u(z  ' )  =  a,  the  width  of  the  current  derivative  pulse 

o  c 

p    '.   Then  from  (17),  a=A(Az')2or 

o 


Az'  =  (a/A)  1/2  (18) 

Using  this  value,  the  maximum  magnetic  radiation  field  for  the 
rising  front  of  the  magnetic  field  pulse  is  easily  evaluated 
from  (13)  under  the  assumption  that  s  and  R  are  slowly 
varying  to  give 


v2  s  ,a^   1/2 

B     =  —  P  '   2   -  (19) 

max  _  2  m     I .. 

cc  R z  x  A 

o 


Values  of  B  for  the  rise  up  to  the  peak  value  given  above 
are  found  by  the  same  process  but  using  appropriately  smaller 
values  of  Az ' .  The  result  is  that  the  integral  (and  therefore 
B)  increases  as  t1/2  after  the  onset  of  the  pulse.  After  the 
maximum  magnetic  field  is  reached,  the  integral  splits  into  two 
parts.  If  the  expression  on  the  right  side  of  (19)  is  called 
1(a),  the  value  of  B  at  later  times  becomes 
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B  =  I (a'  +  a)  -  I (a' )  (20) 

where  a'  is  the  distance  by  which  the  minimum  in  u(z')  is  below 
the  lower  step  of  the  p  '  pulse  in  Fig.  1.   The  first  term  in 
(20)  increases  slowly  with  a',  but  the  second  term  decreases  as 
(a')1/2  leading  to  the  sharp  fall  off  of  the  magnetic  field 
after  the  maximum  as  shown  in  Fig.  2. 

A  complete  current  pulse  may  be  considered  as  a  linear 
rise,  followed  by  a  constant  current,  and  then  a  linear 
decrease.   The  latter  part  gives  rise  to  a  negative  p'(u)  and  a 
reversed  magnetic  field  pulse  so  that  the  magnetic  field  for  a 
complete  current  pulse  has  the  double  peaked  structure  shown  in 
Fig.  2. 


ELECTRIC  RADIATION  FIELD 

In  a  manner  similar  to  the  derivation  of  (13),  the  electric 
radiation  field  may  be  found  from  (2),  (8),  and  (9).   The 
details  are  omitted,  but  the  result  is 

*--f-J(ff--   f)  i    %'<u)dz'         (21> 

0  J  0         0 

The  direction  of  E  may  be  determined  from  the  following 
considerations.   If  R  is  assumed  approximately  constant  and 
denoted  by  R   in  the  reaion  which  contributes  most  stronqly  to 
the  integral,  then 


E--=   E  •  ~   =1(1-  ^cos9)  (22) 

m 


where  I  represents  the  integral  in  (21)  without  the  factor  in 

± 

parenthesis  and  9  is  the  anqle  between  R   and  the  z  axis. 

3  m 

But  the  value  of  R   in  the  region  which  contributes  to  the 

m  3 

integral  is  found  by  evaluating  the  general  expression  (21)  at 

z'  =  z  '.   To  simplify  the  expression,  let  the* observer  be  at 

z  =  o  and  also  assume  that  the  p  '  pulse  is  centered  near 

o 

u  =  o.   Then  R   =  (s2  +  z'2)l/2    may  be  evaluated  using  (14)  to 
m 


give 


("S 


2  x   -1/2 


Rm  =   sfl |  (23) 


10 


If  the  usual  Cerenkov  angle  is  defined  as  cos  9   =  —    R 

c    v     m 

can  be  written  as 


R    =  -r-2-r-  (24) 

m    sin  9 
c 

From  (22 1  it  is  apparent  that  E  is  perpendicular  to  R   when 
y  m 

9=9   and  (24)  shows  that  R   is  the  value  of  R  at  the  Cerenkov 
c  in 

anqle  9  .   Thus,  the  electric  field  from  the  front  of  the  pulse 
c 

(i.e.,  z '  =  z  )  is  transverse  to  R  . 
m  m 

The  situation  is  clarified  in  Fig.  3.   The  charge, 

traveling  from  left  to  right,  emits  a  signal  from  A,  which 

travels  to  the  the  observer  at  0,  traversing  a  distance  R  . 

m 

The  observer  is  at  z  =  o  and  a  distance  s  from  the  path.   The 

field  E  is  oerpendicular  to  R  .   The  signal  was  emitted  from  A 

r  m 

at  an  earlier  time  t'  in  order  to  arrive  at  0  at  the  time  t, 

with  c(t-t')  =  R  .   By  the  time  the  signal  reaches  0,  the 

m     J 

particle  is  at  B,  with  D  =  v(t-t').   Then  R  /D  =  c/v  =  cos9   as 

expected.   D  is  the  path  length  from  A  to  B . 

From  Fig.  3,  one  should  also  note  that  the  field  is 

transverse  to  R  ,  which  points  from  the  earlier  (retarded  time) 

m 

position  of  the  particle,  and  is  radial  relative  to  the  present 

4 

position  of  the  particle.   The  former  condition  holds  for 
typical  dipole  radiation,  while  the  latter  condition  holds  for 
the  Lienard-Wiechart  field  for  a  particle  moving  with  v  <  c. 
The  Lienard-Wiechart  fields  fall  off  as  the  inverse  square  of 
the  distance,  and  do  not  represent  radiation.   In  contrast,  the 
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fields  discussed  here  fall  off  more  slowly  than  R~ 1  and 
represent  radiation;  and  the  total  radiated  power  is  discussed 
in  the  next  section. 


12 


RADIATED  POWER 

The  enerqy  radiated  may  be  found  bv  calculatina  the 
Povntinq  vector  and  inteqratina  over  a  surface.   If  the  surface 
is  a  cylinder  centered  on  the  z  axis,  the  fields  at  a  qiven 
time  have  a  pattern  independent  of  anqle  and  a  z  deDendence  as 

shown  in  the  top  curve  of  ^iq.  2.   The  Poyntinq  vector  is,  of 

-► 
course,  alonq  R  ,  and  the  outward  comDonent  mav  be  inteqrated 

m 

over  the  cvlinder  to  qive  the  total  oower  radiated.   As  a  crude 
estimate  for  the  inteqral,  reolace  the  field  bv  the  oeak  field 
(19)  and  let  the  spatial  width  of  the  Dulse  be  a.   The  radiated 
power  is  then 


R        oo  ( 25 ) 

P(apDrox)  =   —   vl  2sin29  [ 

c2      o       c 


in  cqs  units.   In  the  mks  system,  the  sauare  bracket  is 

replaced  by  2uo/it. 

5  , 6 
In  the  earlier  calculations,    the  fields  and  power  were 

expressed  in  terms  of  Fourier  amplitudes.   If  the  same  current 

Dulse  is  assumed,  P   has  frequencv  components  up  to  the  value 

of  co  such  that  the  wave  lenqth  of  the  radiation  is  ecrual  to  the 

pulse  lenqth.   If  it  is  assumed  that  P   rises  linearlv  up  to 

this  frequency  and  suddenly  drops,  the  total  oower  radiated 

becomes  (in  Mks  units) 
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71 

2 


P   =  7r    v      v  I  2  sin29 


(26) 


Equations  (25)  and  (26)  are  both  rough  estimates  and  the 
point  is  that  the  similarity  of  the  results  is  asserted  to  be 
confirmation  that  the  calculations  in  this  paper  represent  the 
Cerenkov  radiation,  here  expressed  in  terms  of  time  dependence 
of  the  fields. 


14 


DISCUSSION 

In  preceding  sections,  the  time  structure  of  the  electric 
and  maqnetic  radiation  fields  was  developed.   Only  the  far 
fields  were  retained  and  in  the  development  leadinq  to  the  B 
field  (13),  and  the  E  field  (21),  only  the  assumption  of  a 
riqid  charqe  confined  to  a  line  was  introduced.   It  is  seen 
from  these  equations'  that  the  time  derivative  of  the  current  is 
:he  source  function. 

The  simple  model  chosen  to  demonstrate  the  method  of 

3evelooinq  the  time  structures  was  that  of  a  uniform  charqe 

Jistribution  with  uniformlv  varyinq  front  and  rear  sections. 

[•his  model  qives  the  sauare  pulse  charqe  derivative  of  heiqht 

d1   shown  in  Fiq .  1  which  is  easv  to  use  in  evaluatinq  the 
m 

Lnteqral  (13).   Similar  remarks  hold  for  the  power  series 
expansion  of  u(z')  which  is  an  increasinqly  better 
approximation  as  the  time  durinq  which  the  current  is  chanqinq, 
decreases.   Current  variations  other  than  linear  mav  be  readilv 
Incorporated  within  the  framework  aiven.   Also  it  should  be 
loted  that  in  all  cases  the  variation  of  R~2  in  ( 1  "M  which  was 
issumed  constant  in  the  example  will  tend  to  sharpen  both  the 
Leadinq  and  trailinq  edges  of  a  field  pulse. 

In  the  evaluation  of  the  time  structure  of  the  fields,  the 
Deak  field  arose  when  the  inteqral  (13)  had  the  most  widely 
spaced  limits;  a  situation  which  occurs  because  u(z')  has  a 
leqative  slope  for  sufficientlv  neqative  values  of  (z1  -  z)  as 
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shown  at  the  left  side  of  Fig.  1.   In  the  non-Cerenkov  case 

(—  <  1)  this  situation  does  not  arise  since  then  the  slope  of 
c 

the  u(z')  function  always  has  the  same  sign.   In  this  case 
(i.e.  v  <  c)  the  u(z' )  curve  bends  downward  instead  of  upward 
for  large  negative  values  of  (z'  -  z)  and  the  only  contribution 
to  the  integral  (13)  is  from  small  regions  of  z'. 

These  results  show  how  the  time  structure  of  Cerenkov 
radiation  arises  from  the  time  rate  of  change  of  the  charge 
distribution  in  an  electron  bunch.   Present  technology  is  such 
that  this  structure  is  not  observable  in  the  Cerenkov  radiation 
from  S  or  L  band  Linacs  because  of  their  relatively  high 
fundamental  frequency.   However,  induction  accelerators  with 
their  longer  electron  bunch  structure  should  produce  Cerenkov 
signals  in  air  for  energies  greater  than  about  25  Mev,  which 
should  be  observable. 

The  extension  of  this  method  of  calculation  of  the  fields 
for  both  Cerenkov  and  sub-Cerenkov  charge  velocities  is  easily 
made  for  cases  for  which  the  charge  derivative  p'(u)  is  not 
constant.   A  detailed  report  is  under  preparation. 

Finally,  we  note  that  although  the  results  of  other 

4  7 
workers  '   often  have  singularities  in  the  radiated  power  at 

the  Cerenkov  angle,  the  present  results  and  our  previous 

ones  '   show  that  the  radiated  power  is  finite  whether 

calculated  in  the  frequency  or  time  domain.   Also  it  should  be 

noted  that  causality  is  satisfied  because  the  fields  are  zero 

at  times  earlier  than  the  leading  edge  of  the  pulse  shown  in 

Fig.  2. 
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FIGURE  CAPTIONS 

Figure  1.   The  function  u  =  z'  -  vt ' ,  defined  in  the  text  is 

plotted  for  increasing  times  t  \,    1 1>    t 3  at  the  point 

of  observation.   The  corresponding  current  derivative 

profile,  on  the  right,  is  a  function  of  u  only  and 

remains  fixed  in  time.   The  field  signal  pulse  starts 

at  t  ,  and  reaches  a  maximum  at  t3 . 
2 

Figure  2.   The  electric  field  pulse  generated  by  the  beam 

current  profile,  shown  in  the  lower  curve. 

Figure  3.   Geometrical  relations  for  the  Cerenkov  pulse.   The 

source  ( p  ' )  at  A  emits  a  signal  at  an  early  time 
0 

giving  the  E  field  at  the  observer  9;  when  the  field 
reaches  the  observer,  the  particle  is  at  B. 
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